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Everyone is familiar with the climate phenomenon "El Nino” and the havoc it 
causes to weather around the globe. Scientists, using sophisticated computer models and 
observations of temperature, winds, and clouds are trying to forecast when an El Nino 
will occur. Sometimes, if the conditions are right, enhanced precipitation and winds 
move from the eastern Indian Ocean to the central Pacific and may trigger an El Nino 
event. Models have difficulty predicting El Nino in part because they have difficulty 
predicting these propagating systems. Thus, scientists need to describe these systems 
better. This paper offers one mechanism for the generation of these systems. 

One necessary condition is more than normal rain just south of the Equator in the 
eastern Indian Ocean and less than normal rain centered over the Equator, south of India. 
This rain pattern is consistent with strong winds blowing to the east over Indonesia. 
However, this alone is not enough. These winds must be "pumped" towards the central 
Pacific. This is accomplished by switching the precipitation pattern described above on 
and off, essentially creating these bursts of wind every 30 to 60 days. These conditions 
were observed 8 and 6 months before the two strongest El Ninos in the last 23 years, but 
were not observed with weaker El Ninos. Thus, this mechanism is not foolproof for 
predicting El Ninos, but it adds one piece to the bigger puzzle of El Nino. 




Abstract 


A pattern of variability in precipitation and lOOOmb zonal winds for the tropical Indian 
Ocean during 1979 to 1999 (AtmlO mode) is described using EiOFs. The AtmlO mode consists 
of a cross-equatorial gradient of precipitation anomalies and equatorial wind anomalies of 
alternating signs on the Equator. The positive phase is defined as enhanced precipitation to the 
south of the equator, suppressed precipitation to the north, and anomalous westerlies centered on 
the island of Sumatra. In September-October 1981, February-March 1990, and October- 
December 1996 the AtmlO mode was positive and there was a significant 30-60 day variability 
in the gradient of precipitation anomalies. These cases coincided with moderate to heavy 
activity in the Madden-Julian Oscillation (MJO). Links between the AtmlO, MJO, and El Nino 
are discussed. 

1. Introduction 

Recently, progress has been made towards characterizing internal modes of climate 
variability in the Indian Ocean. Webster et al. (1999) and Saji et al. (1999) observed an 
equatorial dipole of sea surface temperature (SST) with links to winds and precipitation, and 
proposed forcing mechanisms that would maintain this ocean-atmosphere interaction. Although 
this mode may lead to better climate predictions in the Indian Ocean region, it does not appear to 
be related to the El Nino / Southern Oscillation (ENSO) and would not improve the skill of such 
forecasts. This study complements Webster (1999) and Saji et al. (1999) by exploring 
atmospheric modes at higher frequencies. 

Ever since Walker i 1924) recognized the Southern Oscillation in trying to forecast Indian 
monsoon rainfall, there. have been efforts to relate the climate of the Indian Ocean sector to 



climate variations elsewhere (Rasmusson and Carpenter 1983, Webster and Yang 1992, Webster 
et al. 1998, Reason et al. 2000). Although much of the previous work has focused on the 
seasonal to interannual time scales, the eastern Indian Ocean is often a starting point for 
submonthly migrations of convection and the Madden-Julian Oscillation (MJO) (Madden and 
Julian 1994, Meehl et al. 1996). 

From 21 years of Global Precipitation Climatology Project (GPCP) data (Huffman et al. 
1997, Pingping Xie personal communication) and NCEP reanalysis winds (Kalnay et al. 1996), 
this study finds a mode of atmospheric variability in the tropical Indian Ocean sector that is 
distinct from the dipole mode proposed by Webster et al. (1999) and Saji et al. (1999). 

2. Results 

Figures la and lb show the first and second EOFs of monthly precipitation (pEOFl. 
pEOF2) for the Indian Ocean domain (50-105E; 30N-30S). The pEOFi, with negative (positive) 
values in the eastern (western) Indian Ocean, is clearly a representation of the dipole mode, 
resembling Figure 4 in Saji et al. (1999). The pEOF2 explains almost 9% of the variance (Fig. 
lb), and is a significant mode of variability according to the criteria of North et al. (1982). 
While pEOFl is essentially a zonal gradient in precipitation anomalies, pEOF2 can be 
characterized as a Southeast-Northwest (SE-NW) gradient about the Equator. 

The first four EOFs of zonal wind anomalies at lOOOmb (u) were calculated. Table 1 
shows the time series correlations between the wind EOFs and the first two precipitation EOFs 
(Figs. la,b). The temporal correlation between pEOF2 and the fourth EOF of u (uEOF4) is the 
only one that reaches the 1 % significance level, using Quenouille’s (1952) method to account for 
the resolution of effective number of degrees of freedom due to persistence. The time series 


plots for these EOFs are shown in Fig. lc. The correlation improves to +0.26 when precipitation 
leads the winds by a month. Overlain on pEOF2 is uEOF4 (Fig. lb). The basin is dominated by 
negative (easterly) anomalies, with positive (westerly) anomalies south of 25S and at the Equator 
from 75E to 105E and near 50E. The pEOF2 and uEOF4 will be defined as the Atmospheric 
Indian Ocean (AtmlO) mode. Seasonal averages did not show a preferred season for the positive 
phase of the AtmlO mode (Fig. lb). However the negative phase (opposite sign of Fig. lb) is 
most likely to occur in boreal winter. This suggests some overlap between the negative phase of 
the AtmlO mode and the Indian Ocean dipole. 

Fig. 2 shows the absolute magnitude of correlations between Nino 3.4 and the dipole and 
AtmlO modes, when the EOFs lead from 16 to 0 months. The Nino 3.4 autocorrelation is plotted 
as a measure of the skill of persistence. As noted by Saji et al. ( 1999) and Webster et al. ( 1999 ). 
the dipole pattern shows small correlation with ENSO. Although the AtmlO mode has equally 
small correlation, these values are better than persistence at a lead time of 9 to 14 months. 

In order to explore the links between precipitation, equatorial wind, and ENSO in the 
eastern equatorial Pacific, the key feature of the AtmlO mode (the SE-NW gradient of 
precipitation anomalies) was reconstructed at a 5-day time step. The daily NCEP wind data was 
determined to be of lesser quality and thus not used. The SE-NW gradient of precipitation is 
assessed by computing the average anomaly within the Northwest box (NW, Fig. lb) subtracted 
from the average anomaly within the Southeast box (SE, Fig. lb). The time series of 
precipitation anomalies (Fig. 3a) shows a substantial amount of high frequency variability. A 
wavelet analysis (provided by C. Torrence and G. Compo and available at the URL 
http://paos.coIorado.edu/research/wavelets) was performed on this time series (Fig. 3b). An 
inspection of the high frequencies (0.0312 to 0.125 years) reveals significant power (+90/7:) in 



nearly all subseasonal bands at the end of 1981 and 1996. Fig. 4a shows the 30-60 day scale- 
average time series (thin line). Seven times during the record the mean 30-60 day power is well 
above the 95% confidence level. However, these peaks coincide with a strongly positive mean 
state of the AtmlO (thick line) only in 1981, 1990, and 1996. Fig. 4b shows a scatterplot of the 
time series in Fig. 4a. In general low-frequency variability is strong when high-frequency 
variability is weak and vice versa, as shown by the envelope of scatter. However, several points 
falling outside this envelope were identified as occurring in September-October 1981, February - 
March 1990, and October-December 1996 (Fig. 4b). For these cases the state of the gradient was 
positive and there was an enhancement of the 30-60 day variability. All were seasons of 
moderate to strong MJO activity and the 1981 and 1996 events led the onset of El Nino 
(Nino 3.4 > 0.5) by 8 and 6 months respectively. 

3. Discussion 

This study has identified consistent modes of variability between precipitation and 
lOOOmb zonal wind in the tropical Indian Ocean region. The AtmlO mode is different from the 
dipole mode in many respects. The AtmlO mode is a purely atmospheric mode; the AtmlO 
mode has a pronounced meridional component; the AtmlO mode is not always tied to the annual 
cycle; and the AtmlO mode seems to lead the El Nino phenomenon by roughly a year. In the 
following discussion a mechanism linking the precipitation and zonal wind patterns is proposed. 
A physical connection to El Nino will also be suggested. 

In the positive phase of the AtmlO mode the center of enhanced rainfall (Fig. lb) is well 
south of the equator. This is consistent with the zonal wind pattern (Fig. lb) which is positive to 
the north of the maximum precipitation anomaly and negative to the south (clockwise flow'). 


Secondly, the larger area of suppressed rainfall (Fig. lb) sits on the equator and is accompanied 
by anomalous westerlies to the East and easterlies to the West, which is consistent with steady 
equatorial dynamics. It appears that the long-term (six months) average of this pattern of 
precipitation anomalies must be positive, and accompanied by a strong 30-60 day variability for 
there to be a subsequent connection to El Nino. We propose that the oscillation in precipitation 
pattern is accompanied by an oscillation in the sea level pressure pattern which serves as a 
"pump" for the surface westerlies over Sumatra. The time scale is consistent with the MJO, 
which has been observed to propagate from the Indian Ocean to the western Pacific (Weickmann 
and Khalsa 1990). It has been shown that winds in the far western Pacific have predictive skill 
for ENSO events, and in particular El Ninos (Clarke and Van Gorder 2001). Therefore, the 
strength and subseasonal variability of the AtmlO mode could be used to predict (up to a year in 
advance) El Nino events that have their very beginnings in the Indian Ocean sector. However, 
further work is needed. Reliable observations of winds and atmospheric pressure, or appropriate 
regional modeling studies are required to advance our understanding of this mode of variability. 
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Figure Captions 

Figure 1 . The first and second EOFs of precipitation (pEOFl, pEOF2) and fourth EOF of 
lOOOmb zonal wind (uEOF4) anomalies in the Indian Ocean sector. A) pEOFl (explains 18.1% 
of the variance). B) Color shading denotes pEOF2 (explains 8.79c of the variance) and contour 
lines denote uEOF4 (explains 7.2% of the variance). Boxes indicate areas used to define the 
SE-NW gradient in precipitation anomalies. C) Time series of pEOF2 (dashed line) and uEOF4 
(solid line). 

Figure 2. Absolute magnitude of time-lag correlations. Solid line is the autocorrelation of Nino 
3.4 from 0 to 16 months. L.ong-dashed line is the correlation between Nino 3.4 and pEOFl when 
pEOFl leads by 0 to 16 months. Short-dashed line is the correlation between Nino 3.4 and 



pEOF2 when pEOF2 leads by 0 to 16 months. Dot-dashed line is the correlation between Nino 
3.4 and uEOF4 when uEOF4 leads by 0-16 months. 

Figure 3. Analysis of 5-day time series of the SE-NW gradient of precipitation anomalies (see 
boxes in Fig. lb). A) Raw time series of the SE-NW gradient. B) Morlet wavelet power 
spectrum of the raw time series in (A). The contour levels are chosen so that 75%, 50%, 25%, 
and 5% of the wavelet power is above each level, respectively. Cross-hatched region is the cone 
of influence, where zero padding has reduced the variance. Black contour is the 10% 
significance level, using the global wavelet spectrum as the background. 

Figure 4. Relationship between the power of 30-60 day oscillations and long-term running 
average of the SE-NW gradient of precipitation anomalies. A) Thin line denotes time series of 
the 30-60 day scale-average of the wavelet power spectrum in Fig. 4b. Thick line denotes time 
series of the six month running mean of the SE-NW gradient of precipitation anomalies (Fig. 4a). 
B) Scatterplot of the two time series in (A), showing precipitation anomalies as a function of 30- 
60 day power. Certain data points are identified with open circles, squares, and diamonds 

Table 1. Correlations between precipitation and u EOFs 

uEOFl uEOF2 _ uEOF3 uEOF4 
pEOFl -0.09 +0.19 4403 1405 
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a) SE-NW gradient of anomalies (pentad) 
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30-60 Day Scale-average Time Series (Fig. 3b) And 
Six Month Running Mean of Pentad Time Series (Fig. 3a) 
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Scatterplot of Time Series in (A) 
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